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We show how the formulation of the matrix models as conformal field theories on a Riemann sur- 
faces can be used to compute the genus expansion of the observables. Here we consider the simplest 
example of the hermitian matrix model, where the classical solution is described by a hyperelliptic 
Riemann surface. To each branch point of the Riemann surface we associate an operator which rep- 
resents a twist field dressed by the modes of the twisted boson. The partition function of the matrix 
model is computed as a correlation function of such dressed twist fields. The perturbative construction 
of the dressing operators yields a set of Feynman rules for the genus expansion, which involve ver- 
tices, propagators and tadpoles. The vertices are universal, the propagators and the tadpoles depend 
on the the Riemann surface. As a demonstration we evaluate the genus-two free energy using the 
Feynman rules. 



1 Associate member of the Institute for Nuclear Research and Nuclear Energy, Bulgarian Academy of Sciences, 72 
Tsarigradsko Chaussee, 1 784 Sofia, Bulgaria 



1 Introduction 



In this paper, which is a continuation of [1], we formulate the 1/N expansion of U(N) invariant 
matrix integrals as the quasiclassical expansion in conformal field theories on Riemann surfaces. Our 
basic example will be the hermitian matrix model defined by the partition function 



Z N [V] = J dMe- trV W 



N 



C 



Y[dzie- V M Hizi-zj) 2 . (1.1) 

i=l i<j 



Here Zi, . . . , zn are the eigenvalues of the matrix variable M and 

oo 

V(z) = -J2 t nZ n (1.2) 

n=0 

is a confining potential. For the hermitian matrix integral the contour of integration goes along the 
real axis. Then the integral is convergent if the potential becomes infinitely large when z — > ±oo. If 
this is not the case, the integral can be made convergent by deforming the contour of integration. 
The 1/N expansion of the free energy J^fV] = lnZAr[y] is of the form 

oo 

T N [V] = N 2 - 2g F {9) [V/N] + nonperturbative terms, (1.3) 

3=0 

where the term counts the number of genus g fat Feynman graphs. 2 The leading, genus zero, 
term can be obtained by evaluating the saddle point eigenvalue distribution [2]. In this approximation 
the eigenvalues can be considered as a continuous charged liquid defined by the spectral density 

1 - 

p(z)= lim t; (1-4) 

7V->oo N t— 1 
i=l 

The collective field theory for the spectral density ( 1 .4) is ill defined at small distances and cannot 
be used to compute the higher terms of the quasiclassical expansion. This can be done by solving the 
Ward identities for the integral (1.1), known also as loop equations. An improved iterative scheme 
for calculating the higher genus observables, known as method of moments, was set up by Amb- 
jorn, Chekhov, Kristjansen and Makeenko [3]. More recently, Eynard [4] formulated the recursion 
procedure for solving the loop equations in an elegant graphical scheme. 

The form of the loop equations suggests that there is a local field theory hidden behind. Indeed, 
the loop equations can be formulated as Virasoro constraints for the bosonic field 

N 1 
0(z) = V2Y,log(z-z i )-—V(z). (1.5) 

i=i v 2 

Therefore the theory in question must be a two-dimensional chiral CFT An important feature of this 
CFT is that the bosonic field develops a large classical expectation value <fi c {z)- The classical current 



"Strictly speaking, this is true only in the case of one-cut solutions. 
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Jc = d z 4> c satisfies a quadratic equation, the classical Virasoro constraint, whose solution defines for 
a general potential a hyperelliptic Riemann surface with branch points a± , . . . , a2p placed on the real 
axis. The set of branch cuts along the segments [01,02], ... , [a2p-i, <J2p] coincide with the support of 
the classical spectral density. Each branch cut is associated with a local minimum of the potential and 
is characterized by the filling number Nj, equal to the number of the eigenvalues trapped around this 
minimum. 

In the chiral CFT the branch points of the Riemann surface can be thought of as the result of 
applying twist operators to the sZ(2) invariant vacuum. The twist operators are however not conformal 
invariant. In particular the translations change the positions of the branch points. In [1] we claimed 
that in order to satisfy the conformal Ward identity the twist operators should be dressed by the modes 
of the twisted boson. The dressed twist operators were called in [1] star operators because of the 
analogy with the star operators introduced by G. Moore in [5]. We should stress here that the star 
operators in [1] and in [5] are different objects. 

The problem of computing the genus expansion of the of the observables is thus reduced to the 
perturbative construction of the operators dressing the twist fields. Such dressing operator represents 
a formal series expansion in the modes of the twisted boson. The coefficients of the expansion are 
related to the correlation functions in the Kontsevich model. In this sense the proposal made in [1] 
allows to decompose the hermitian matrix model with an arbitrary potential into 2p Kontsevich models 
coupled through the modes of a gaussian field on the hyperelliptic surface. Similar decomposition 
formulas for the hermitian matrix model were suggested in [6, 7, 8]. 

The aim of this paper is to present the details of the construction of the dressed twist operators and 
to set up the Feynman diagram technique for calculating the free energy and the expectation values in 
all orders on 1/JV. 

The paper is organized as follows. In Section 2 we remind the Fock space representation of the 
integral (1.1) in terms of a chiral bosonic field <p(z) with Liouville-like interaction given in [1] and 
derive the loop equations from the conformal Ward identity. This Section helps to make the presen- 
tation self-consistent and can be skipped by the reader who is familiar with the loop equations. In 3 
we give the representation of the leading and the subleading orders of the free energy in terms of the 
correlation function of 2p twist operators. In Section 4 we give the explicit prescription for dressing 
the twist operators and set up the Feynman rules for computing all orders of the genus expansion. This 
Feynman diagram technique is different than the diagram technique in [4]. In Section 5 we compute, 
using the Feynman rules, the genus two free energy and compare with the result of [3]. 

2 The matrix model as a chiral CFT 

2.1 From Coulomb gas to compactified chiral boson 

The Coulomb gas integral (1.1) can be represented as a Fock space expectation value in a theory of a 
Liouville-like chiral CFT with the wrong sign of the kinetic term. Introduce the holomorphic scalar 
field 4>{z) with mode expansion at z = 00 

4>(z) = q + J lnz-^— z~ n , [J n ,Jm]=n5 n+m ,o, [J ,q] = l, (2.1) 

and left and right Fock vacua defined by 

(0|J = <?|0) =0, (0|J_„ = J n \0) =0 (n>l). (2.2) 
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The field has operator product expansion 

<j)(z) 4>{z') ~ ln(z - z') (2.3) 
and give a representation of the su{2) current algebra: 

J = ^ dtf>(z), J± =: e ± ^ : . (2.4) 
Then the partition function (1.1) is given by the scalar product 

Z N = (N\e J[v] e Q +|0), (2.5) 



where 



/dz 
—V(z)J(z) = ^t n J n , (2.6) 



OO 
OO 



n>0 

Q+ = / dz,J + (z) (2.7) 

and (N\ is the charged left vacuum state 

{N\ = (Ole^*? . (2.8) 

The operator generates screening charges, the operator e J ^ produces the measure e~ v ^ for 
each charge, and the left vacuum (N\ projects onto the sector with exactly N screening charges. 
Similar Fock space representations of the eigenvalue integral have been proposed in [9, 10, 11]. 
The currents J, J+ and J_ are invariant under the discrete translations of of the field 4>: 

4>^4> + my/2. (2.9) 

Therefore the the field cf> can be compactified at the self-dual radius R s .d. = l/\/2- Furthermore, the 
transformation 

$-*■-</> (2.10) 

is an automorphism of the current algebra. The geometrical meaning of the symmetry (2.10) will 
become clear when we consider the quasiclassical limit of the bosonic field. 

The correlation functions of in the matrix model are obtained through the identification 

4>(z) = V2tr log(z-M) - -j=V{z) (2.11) 

v2 

or, in terms of the su{2) currents, 

J(z) = tr — ^ - \V'(z) , J±(z) = e^[det(z - M)] ±2 . (2.12) 
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2.2 Conformal Ward identity 

To prove that the theory is conformal invariant we have to demonstrate that the energy-momentum 
tensor 

T(z) = \ : d<f,(z)d<f>(z) : = ^ L n z~ n - 2 (2.13) 
commutes with the screening operator Q + . Indeed, for any n > —1 we have 

/OO /* OO ^3 

dz[L n ,J + (z)} = dz — [z n+l J + {z)). (2.14) 
'OO >/ — OO 

Since our potential diverges at infinity, the boundary terms vanish and the result is zero for all n > — 1. 
As a consequence, the expectation value 

(T{z)} = (N\e J[v] T{z)e Q +\0) (2.15) 

is regular for z ^ oo. This condition can be written as a contour integral which projects to the positive 
part of the Laurent expansion of T(x): 

dz' /T(z)-T(z')\ 

™( U ) = a <2 - 16) 

The conformal Ward identity (2.16) is translated into a set of differential Virasoro constraints on 
the partition function using the representation of the gaussian field as a differential operator acting on 
the partition function, 

*M - **> = T2 £ inZ " + ^ ,nz w a + ^ < 2 - 17) 

n>0 n>0 

The Virasoro constraints read 

L n -Z N = (n>-l) (2.18) 

where 

^ d =E^7^ + I><^' ^-Z N = NZ N . (2.19) 
f^dt k dt n _ k f-^ dt n+k dt 

k=0 fe=0 



3 The quasiclassical limit 

3.1 The classical solution as a hyperelliptic curve 

Applied to the genus expansion (1.3) of the free energy, the Virasoro constraints (2.18) generate 
an infinite set of equations for the correlation functions of the current J(z), which can be solved 
order by order in 1/N. The lowest equation is the classical Virasoro condition, which determines the 
expectation value J c (z) of the current in the large N limit. In our normalization J c is of order of N, 
just as the confining potential V. 
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Figure 1 : The A and B canonical cycles for a genus 2 spectral curve (p = 3) 

The classical Virasoro constraint states that T c = J 2 is an entire function of z. The most general 
solution is 

2p 

J c (z) = -M(z) y(z), y 2 = ]\{z - a 3 ). (3.1) 

where M(z) is an entire function of z. Assuming that a\ < a2 < ... < a2 P , the meromorphic function 
J c (z) is discontinuous along the intervals [a,2k-i, a>2k] an d its discontinuity is related to the normalized 
classical spectral density by 

J c (x - jO) - J c (x + tO) 

Our aim is to construct a CFT associated with the classical solution J c . It is advantageous to think 
of the gaussian field as defined not on the complex plane cut along the intervals [a2fc-ij Q2fc]> t> ut on 
the Riemann surface representing a two-fold branched cover of the complex plane, the two sheets of 
which are sewed along the cuts [02^-1, 023]- m this way we trade the boundary condition along the 
cuts for the monodromy relations (cj) — > —<ft) when one moves around the branch points. 

The hyperelliptic Riemann surface is characterized by a set of moduli, associated with the canon- 
ical A and B cycles. The cycle Aj encircles the cut [02^-1 j 0-2 j] and the cycle Bj encircles the points 
a2j, 02^-3, passing through the j-th and the p-th cuts (Fig. 1), so that 

A k oB 3 = 5 Kj (j = l,...,p-l). (3.3) 

The classical current is determined completely by the potential V(z) through the asymptotics 

J c (z) = -\V\z)+Nz- 1 + ... (3.4) 

and by the charges N\ + • • • + N p = N associated with the ^4-cycles, 

dz 

— J c {z) = N 3 , j = l,...,p. (3.5) 



From (3.5) it follows that the derivatives 



dJ c( z ) a 

Uj = -^-dz (3.6) 



5 



form a basis of holomorphic differentials of first kind, associated with the cycles A^, 




(3.7) 



The integrals of ojj along the S-cycles give the period matrix of the hyperelliptic curve 

r kl = ±-J^. (3.8) 

In the following we will consider the filling numbers Nj as fixed external parameters. The cor- 
responding free energy is that of a metastable state, which becomes stable up to exponentially small 
effects in the large N limit. We will understand the 1/N expansion (1.3) of the free energy in this 
setting. Alternatively, one can introduce chemical potentials Tj for the charge Nj and evaluate the 
partition function for fixed Fj . 

If one is interested in the quasiclassical evaluation of the original partition function (1.1), one 
should perform the sum over all possible Nj. This problem has been considered and solved in [12]. 
After performing the sum, the logarithm of the partition function does not have the 1/N expansion 
(1.3). Since the numbers Nj are the discontinuities of the bosonic field along the A-cycles, the sum 
over all Nj means that the bosonic field is effectively compactified at the selfdual radius. 



3.2 The branch points as primary conformal fields 

Our goal is to construct a CFT associated with the classical solution J c . We will think of the hyper- 
elliptic Riemann surface as the complex plane with conformal operators, the twist operators with 
dimension 1/16, associated with the branch points. This point of view was first advocated by Alexei 
Zamolodchikov in [13]. Later some of the findings of [13] were obtained independently by Dixon at 
al [14] and further developed by a series of brilliant papers by V. Knizhnik [15, 16, 17]. 

Let a be one of the branch points of the classical solution J c . In the vicinity of a the current has 
mode expansion 

J(z) = J c (z)+ Y, Jr{z-aY r ~ l (3.9) 

with the following algebra, 

[J r , J s ] = \r5 r+sfl . (3.10) 

The Ramond vacuum associated with this branch point is defined as the highest weight vector of 
the representation of this algebra. The corresponding quantum field is the twist operator a (a) 

Jx/ 2 o-(a) = J 3 /20-(a) = J 5 / 2 o-(o) = • • • = 0. (3.11) 

The HUbert space associated with a (a) generated by multiple action on a (a) with the negative modes 
J_l/2; J-3/2i • • • of the current J. 

The product of two twist operators a{a\)a{a2) is a single-valued operator with respect to the 
current J(z). This state can be decomposed into charged eigenstates (3.5). Let us denote a state with 
given charge N by [cr(ai)a(a2)]N- Our strategy in the following is to simulate the cuts [d2j-i, a-2j] of 
the Riemann surface by multiplying the right vacuum with the operators [a(a,2j-i)(r(a2j)]N ■ an d give 
a Fock-space representation of the partition function of the matrix model similar to (2.5). In the new 



6 



representation the operator creating the N charges is replaced by a product of p pairs of twist 
operators. 

The correlation function of p pairs of twist operators was calculated by Al. Zamolodchikov [13] 
in the case when the current has no expectation value and the total charge is zero: 

(0| n^-iVMklO) = Z twist ( fll ,...,a 2p ) j^lk=i^ N i N * . (3 . 12 ) 
The meromorphic function Z tw i st (ai, . . . , a2 P ) is given by 

2p 

Ztrt* = n {fli - a k y 1/8 [detK]- 1 / 2 , (3.13) 

j<k 

where the matrix K is denned by 

r z i~ l dz 

K v= ~ TT (i,j' = l,-..,p-l)- (3-14) 

It is straightforward to generalize this formula to the case of non-zero total charge and non-vanishing 
expectation value of the current. We should simply insert the operator exp J[V] as in (2.5). This gives 
a representation of the partition function in the quasiclassical limit as the scalar product 

quasicl = ^ e J [V] f[ [a ( a2j _ l)aM]N .\ y (3.15) 

The normalized expectation value of the current evaluated with respect to (3.15) obviously coin- 
cides with the classical solution J c : 

( J{*) > q uasicl = 75 dj>{z) ■ log^ UaSid = MZ) • (3.16) 

Therefore the r.h.s. of (3.15) satisfies the classical Virasoro constraints and reproduces correctly the 
leading order of the free energy. The quantum Virasoro constraints are however not respected. The 
twist operator a (a) depends on the position of the branch point a and does not satisfy the lowest 
Virasoro condition L_i. It is also not invariant with respect to dilatations generated by Lq. 

4 The genus expansion 
4.1 Dressed twist operators 

We would like to modify (3.15) so that it reproduces all orders of the l/N expansion of the free energy. 
We have to look for operators which create conformally invariant states near the branch points. Such 
operators can be constructed from the modes of the twisted bosonic field near the branch point by 
requiring that the singular terms with their OPE with the energy-momentum tensor vanish. 

In [1] the author proposed 3 that twist operator a (a) can be made conformal invariant by mul- 
tiplying with an appropriate dressing operator e w( - a \ which is made from the modes of the twisted 
field: 

a(a) -»• 5(a) = e^ a) a{a). (4.1) 
3 This proposal was reproduced in more details in sect. 4.3 of [18] following an unpublished extended version of [1]. 
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Assuming that the dressed twist operator is a well denned operator in the Hilbert space associated 
with a(a), the dressing exponent w(a) can be expanded as a formal series in the creation operators 
J- r (r > 0) defined by the expansion (3.9) near the point a: 

™( a ) = J2^ w ri ...r n (a) J- ri -J-r n - (4-2) 

n>0 ' ri,...,r n 

The coefficients of the expansion depend on the classical current J c and are determined by the re- 
quirement of conformal invariance 

L n S(a) = (n > -1) (4.3) 

where the Virasoro generators L n are defined by the expansion at the point z = a of the energy- 
momentum tensor 

= Y J L n (a)(z-a)- n - 2 . (4.4) 

n 

One finds for the Virasoro operators with n > — 1 

L n = ■ (Jr + Jr) (Js + J!) ■ +^n,0, (4-5) 
r+s=n 

where J'f are the modes in the expansion of the classical current near the point z = a: 

J c (z)= ]T Jl r (z-a) r - 1 . (4.6) 

r>3/2 

By convention J c _ r = if r < 1/2. 

Once we found conformal invariant operators that create the branch points, it is clear how to repair 
(3.15) so that it holds for all orders in 1/N 2 . It is still possible to assign to the pair of dressed twist 
operators associated with the endpoints of the cut [02^-1, d2j) a definite charge Nj. We denote this 
state by [S(a 2 j-i)S(a 2j )] Nj . 

Our claim is that the partition function of the matrix model is equal, up to non-perturbative terms, 
to the expectation value 

v 

Z N = (N\ e J [ y l ]\[S{a2j-i)S{a 2j )] Nj \ti), (4.7) 
i=i 

where J(z) is the current of the Z^-twisted gaussian field defined on the Riemann surface of the 
classical solution. Indeed, this expression satisfies the conformal Ward identity in all orders of 1 /N 
since by construction the energy-momentum tensor (4.4) commutes with the dressed twist operators. 
Furthermore, the leading order expectation value of the current J c {z) satisfies the asymptotics at 
infinity (3.4) and the conditions (3.5). 

Remark: The Fock space realization (4.7) of the partition function resembles the QFT represen- 
tation of the r-function for isomonodromic deformations obtained by T. Miwa [19] and revisited by 
G. Moore [5]. In our case there is an irregular singularity at infinity and 2p regular singularities at the 
branch points. For matrix ensembles with hard edges (infinite wall potential), it was convincingly ar- 
gued in [20] that the Fock space representation (2.5) can be transformed into to a correlation function 



T{z) 



lim 



J(z)J(z') 



2(z 



y'\2 
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of Moore's star operators. However, even in this simplest case, the perturbative, or l/N, expansion 
of the star operators is ill defined. The ambiguity gets even worse for smooth potentials. The star 
operator is defined in [5] by the exponential of a contour integral starting at the branch point, but in 
the case of a smooth potential the position of the branch point should be adjusted at each order in 
l/N. The dressed twist operators (4.1) generate the complete perturbative expansion, while the star 
operators of [5] capture the leading non-perturbative behavior. 



4.2 Fock space representation of the partition function 

The l/N expansion can be obtained by considering the dressing operators as a perturbation and ex- 
pand them in the negative modes of the current. Let us denote by the modes of the current J 
associated with the expansion around the branch point aj . Then the total dressing operator, which we 
denote by Cl, is given by the formal series 

a= n ^ (aj) < E «&U J-i-^i- (4-8) 

j=l,...,2p n>0 ' ri,...,r„ 

The partition function Zn is equal, up to non-perturbative collections, to the normalized expectation 
value of the dressing operator with respect to the left and right states 

(left) d =l f (N\ e J [ y l, (right) d ^ J[\a{a 2j ^)a{a 2j )] Nj |0>. (4.9) 

i=i 

In order to perform the expansion we should evaluate the expectation value of any product of 
negative modes Jj? . Since J(z) is the current of a gaussian field, it is sufficient to calculate the 
expectation value of a pair of such modes: 

\a-a-] , Fo-1 fa-1, def (left I I right) 

GrV = (</_;. J™) = - — /.v . (4-io) 

(left (right) 

The matrix G^, can be computed knowing the two-point function ( J(z)J(z') ), which is the 
unique function defined globally on the Riemann surface and having a double pole at z = z' with 
residue 1/2. From the definition (4.10) and the mode expansion (3.9) it follows that 

J 2iri J 2m (z - Oi) r (zf - atf' ' K ' ' 



Once we know the matrix G^s ^ and the coefficients w L r 1 1 }.. rn , we can compute the l/N expansion 
to any order just by expanding the dressing operators and performing Wick contractions. 

This prescription can be packed in a concise formula in the following way. Introduce, together 
with the right Fock vacuum associated with the 2p twist operators, a left twisted Fock vacuum, which 
annihilates the negative modes of the current. The two vacuum states, which we denote by (0 tw | and 
|0tw), are defined by 

(0 tw |jl^ ] =0, J^ ] |0 tw ) = (r > 1/2; j = l,...,2p). (4.12) 

Then the state (right) can be identified with |0 tw ) and the state ( left | is obtained from (0 tw | by acting 
with the gaussian operator associated with the matrix (4.10). As a result we obtain the following Fock 
space representation of the expectation value (4.7), 

Z N = (left | £l (right) = (left(right) (0 tw | e 2jdj Cl |0 tw ) , (4.13) 
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where Cl is defined by (4.8) and 

Jdj = E E — Cte* aj] 4 ai] Jl aj] . (4.14) 

*,i=l r,s>l/2 

4.3 The dressing operator and the Kontsevich integral 

In order to make use of this formula we need the explicit expressions for the coefficients of the series 
(4.2), which can be obtained by demanding that the dressing operator Cl = e w solve the Virasoro 
constraints generated by the operators (4.5). 

A detailed analysis of the solution of these Virasoro constraints can be found in [21]. The solution 
is given by the Kontsevich matrix integral [22], also known as matrix Airy function. To make the 
connection with the notations of [21] we represent the modes J r as 

J-n-i/2 = ~\t n , Jn+i/2 = -{n + \)d n (n > 0) (4.15) 
where we denoted d n = d/dt n . We also relabel the modes of the classical current as 5 

J c _ n _i -+-\nn (n>0). (4.16) 
Then the dressing operator (4.2) is represented by a function Q(to, h, £2, ■ ■ ■ ) which satisfies 

L n n = (n>-l) (4.17) 

with 

L n = (k + ±)(t m + n m )d k + Yl {k + \){m + \)d k 8 m (4.18) 

k—m=n k+m=n—l 

+ 7^1,0 + -^ (n>-l). (4.19) 
4 lb 

The solution depends on the moments fj, n of the classical current. It is sufficient to have the 
solution Oo f° r the simplest nontrivial classical background fj, n = fit 5 n ^. Then the general solution 
is obtained simply by a shift t n — > t n — \i n , n > 2: 

n = exp(-J2vnd n ) fio- (4-20) 

n>2 

The function Uq is given by a formal expansion in io,ii, • • • and 1/fJ-i ~ 1/N: 



g>0 n>0 fci,...,fcn>0 



The coefficients k are the genus g correlation functions in the Kontsevich model. They are 

proportional to the intersection numbers in the moduli space of Riemann surfaces of genus g with n 
punctures: 



,(9) 

i=i 



= (-1)" H^kj - 1)!! {r kl ...T kn ) g . (4.22) 



4 The times t n here are shifted with respect to the times in [21] as (i,i)here = [(2n — 1)!!] 1 (t n — <5n,i)ihcre- 
5 The moments J„ of [21] are related to fi n by 7i = 1 — /ii and I n = — (2n — 1)!! fj, n for n > 2. 
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The intersection numbers (r^ . . . Tfc n ) 9 are positive rationals. They are nonzero only if the indices 
fei , . . . , k n obey the selection rule 

n 

kj = 3( 5 - 1) + n. (4.23) 

The genus zero intersection numbers are the multinomial coefficients 

(Tmi-TmJn = — + ,' ' ' + 7^ ' , mi + ... + m n = U - 3. (4.24) 

m\\...m n \ 

The first several intersection numbers of genus g = 1, 2 are [21] 

(n — 1)! 1 1 

( r Di = — — » ( r o r n+i)! = (Ton^! = Tg, 

7 29 1 

(T9), = , (T2T3), = , (t 4 )„ = . (4.25) 

\ 2/2 240 w 72 5760 x 12 1152 

An efficient procedure for evaluating (r/ Cl . . . t/ Cii ) 9 was proposed in [23]. 
4.4 Feynman rules 

Now we have all the ingredients needed to construct the 1/N expansion of the free energy associated 
with the classical solution with p cuts. Our starting point is the Fock space representation (4.13) of 
the all genus partition function. The expression for the partition function depends on the moments of 
the classical current at the branch points 

and the matrix gm]n^ associated with the two-point function of the current on the Riemann surface 
and defined by 

r K%i = a[*l[*L , 4 21) 

J 2ni J 2iri (z - ai)™+ l l 2 {z' - aj ) n + l l 2 ' K ' 

With each branch point we associate a set of coordinates and use the representation (4.15) 

of the modes of the bosonic current in terms of the Heisenberg algebra generated by tn^ and dn- 
Together with the explicit solution (4.21) for the dressing operators associated with the branch points, 
this leads to the following expression for the genus expansion of the free energy: 



^ ( o )+ ^) = g i7r ££*=i T jk N J N k +Y,„ t„J„ JJ L^Y^Z^,...^)- (4.28) 



2p 

i=i 

2p \ ( 2p 

e E s > 2 iV-^) = exp [ ± £ £ G [ ^ ] dtW exp £ 

i,j=lm,n>0 J \ j=l n>0 



HEE E K" J ) <U J ^r^ • ^ 

K g>0 n>0 fci,...,fc„>0 / t (.) =Q 
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(g) G ^ [ft.*] 



Z' Z II m 



<W(z)> <W(z)W(z')> GJ.z) 



Figure 2: The Feynman rules for genus expansion of the hermitian matrix model. The vertices Wmi,...,m„ 
are universal, the propagators Gm,n depend only on the moduli of the spectral curve, and the tadpoles ' 
are determined by the expansion of the classical solution at the branch points. The vertices are connected to 
each other by propagators and the tadpoles are connected directly to the vertices. The Feynman graphs for the 
correlation functions of the resolvent contain three more elements: the disk, the cylinder and the external legs 
depicted in the second line. 

Using this formula one can evaluate the genus g free energy as a sum of a finite number of Feyn- 
man graphs with vertices kn given by (4.22), tadpoles fin 3 ^ and propagators Gm%^- Note that 
while the propagators and the tadpoles depend on the classical solution, the vertices are universal. 
From (4.25) we find the first several vertices (4.22): 

1 (i) 1 (i) _5 

w n n n — J-i — TTT? iw i 1 — ^TTi w n 9 — ~j w q 12 — Ti w 3 — o' 

4 ' ' o 



,(0) 

0,0,0 — 


"I, 




1 

"24' 


«#= 


1 (1) 

24 °' 2 


,« - 

0,0,2,2 — 


3 
2 


(2) 
W 2,2,2 — 


63 
~80' 


(2) 
W 2,3 = 


29 (2) 
128' 



^0,0,2,2 — 2 W 2,2,2 — W 2,3 ~ W 4 ' ~ ^0, 0,0, 0,2 — 3. (4.30) 

The correlators of the resolvent of the random matrix are obtained from the correlations of the 
current J by the identification (2.12). To evaluate the correlation functions of the current we use its 
representation as differential operator 

2p 

j{z) = j c (z) + ^2Y1 G " j] (*) d ™ jl ■ < 4 - 31) 

j=l n>0 

We therefore extend the set of Feynman rules by adding the external lines which represent the func- 
tions 

The Feynman rues for evaluating the l/N expansion are given in Fig. 2. The genus g contribution 
for any observable is equal to the sum of all genus g Feynman diagrams. The genus of a Feynman 
diagram is equal to 2 — 2h — n, where h is the number of handles, including the handles made of 
propagators, and n is the number of the external lines. 
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5 Example: the single cut solution 

In this section we consider in details the case of one cut (p = 1), where the Riemann surface of the 
classical solution is a sphere. The explicit expression for the classical current is 



Jc{z) 



dz y(z) V'(z') - V'(z) 
j 2iri y(z') z — z' 



l(z) = y/(z- ai)(z - a 2 ). 



(5.1) 



Expanding at z = oo and using the asymptotics J c {z) ~ — \V (z)+N / z+ one finds the conditions 



dz z k V'(z 
A\ 2vri y(z) 



-2N5, 



(A; = 0,1) 



(5.2) 



which determine the positions of the two branch points. 

The two-point function of the ^-twisted current J(z) on the Riemann surface of the classical 
solution is equal to the 4-point function of two currents and two twist operators: 



{J(z)J(z')) ee {0\J(z)J(z')a(a)a(b)\0) 



(z-a 1 )(z / -a 2 ) 
{z-a,2)(z'-ai) 



+ 



(z'-a 1 )(z-a 2 ) 
(«'— 02)(«— ai) 



4(z - z') 



i\2 



(5.3) 



The coefficients G^ ak ' are obtained by expanding {J(z)J(z')) near the points a\ and a 2 . Assuming 
that a,2 > oi, we write 



U km 



dz 



dz 



2-ni L . 2iri 



3 - (z-<ti) k 2(z'- aj y 



(J(z)J(z>)). 



(5.4) 



We find that G^' are of the form 



G 



Oi,Ofc] 
k,m 



Ak+m+1 9km ' 



where g^f"^ are rational numbers with the symmetry g^ ai 
first several coefficients g^lf"^ we 



d=\a\ — a 2 |, 

V 

9km 



(5.5) 



02,021 j r a l> a 2l [08)01 ] T.I 

md 9 km = 9 km • The 



r («i ,«i ) i _ r 1 02,02 n 



(12,02)1 



r (01,02)1 _ r (02,01)-, 



-1/2 


3/8 


-5/16 


3/8 


-3/8 


45/128 


-5/16 


45/128 


-45/128 


-1 


3/2 


-15/8 


3/2 


-21/4 


165/16 


-15/8 


165/16 


-1745/64 



V 



/ 



(5.6) 



As an illustration of the Feynman diagram technique we will evaluate the free energy up to genus 
two. We denote the two branch points by a' = ax, a = a 2 , and the moments of the classical solution 
associated with them by \i n = fjffl , p! n = fi^ ' . 

The genus-one term is 



= 7 In /ii In u'i m<i (d = a — a!). 

24 p 24 P1 8 v ' 



(5.7) 
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Figure 4: The Feynman graphs contributing to the genus two free energy in the case of a single cut. 
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The first two terms come from the dressing of the twist operators and the last term is the logarithm of 
the correlation function (0\a(a')a(a)\0). 

The term J 7 ^ 2 ' in the genus expansion of the free energy is a sum of the contributions of all 
possible genus two Feynman diagrams composed by the vertices and the propagators shown in Fig. 3. 
The relevant diagrams are depicted in Fig. 4. The result is 6 

(2) If 21 ii\ 29 ^2M3 35 fu 5 /i 3 49 \i\ 105 \i 2 175 1 \ 

~i4 \ 160 + 128 nl ~384/I7 + 32^d~256^fd~ 512 d 2 ^ ~ 1024 d? ) 

W l} MiMi V 128 Mid 2 128 ^d 2 256 d 3 J ( J 

The genus two free energy was first computed by Ambjorn et al in [3]. To compare with the result 
of [3] we have to express the moments fi n and fi' n of the classical field in terms of the ACKM moments 
M n and J n , denned as 

2 f dz Mz) 1 (2k - 1)!! 

71 ~ ~N J 2tH (z - a)*+V*{z - a') 1 / 2 ~ N ^ k\(-2d) k ' 

j _ 2 / <** W , r iiy (^-l)» , (59) 

iV jr 27ri (z — a) 1 / 2 (z — a') ,l+1 / 2 1 j iV ^ k\{~2d) k M ™- fe ' k j 

Substituting (5.9) in the answer found in [3] (note that some of the coefficients are corrected in [24]) 

, M 181 181 181J 2 I8IM2 3J 2 



480 J 2 d 4 480M 1 2 d 4 480 J 3 d 3 480M 1 3 d 3 64 J 2 Mid 3 
f 2 HJf 11M| 43M 3 



64Af 2 Jid 3 40,/fd 2 40M^ 2 192Af 3 d 2 
43J 3 J 2 M 2 5 21J 2 3 29J 2 J 3 35J 4 



192J 3 d 2 64J 2 M 2 d 2 16JiMid 4 160Jfd 128Jfd 384J 3 d 
21Af 2 3 29M 2 M 3 35M 4 



160M x 5 d 128M^d 384Af 3 d 
one reproduces the expression (5.8). 



(5.11) 
(5.12) 
(5.13) 



6 Discussion 

The CFT formalism was developed before for the continuum limit of a class of matrix models which 
reduce to Coulomb gas integrals similar to (1.1), such as the SOS and ADE matrix models [25]. 
The classical solutions considered in [25] correspond to non-hyperelliptic Riemann surfaces with one 
higher (or even infinite) order branch point at infinity and one simple branch point on the first sheet. 
The CFT formulation of these models provides a rigorous derivation of the Feynman rules for the 
genus expansion obtained previously in [26, 27]. 

In this paper we developed the CFT formalism for any classical solution of the hermitian one- 
matrix model. We do not see conceptual difficulties to generalize it to any classical solution of the 
above mentioned models, including the 0(n) model. 

The Feynman rules that follow from the CFT representation look similar to the diagram technique 
developed in a larger class of matrix models by Eynard and collaborators [4, 28, 29, 30, 31]. If there 

6 The author thanks A. Alexandrov for pointing out a missing term in the unpublished extended version of [1]. 
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exists an exact correspondence between the two formalisms, then the CFT formalism can be possibly 
extended also for matrix models which do not reduce to Coulomb gas integrals. It is likely that the 
recipe for transforming the Eynard formalism into an effective field theory proposed by Flume et al. 
[32] is the way to obtain this correspondence. 
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